Ferromagnetism and phase separation in one-dimensional 
d— p and periodic Anderson models 
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Using the Density Matrix Renormalization Group, we study metallic ferromagnetism in a one- 
dimensional copper-oxide model which contains one oxygen p-orbital and one copper d-orbital. The 
parameters for the d-p model can be chosen so that it is similar to the one-dimensional periodic 
Anderson model. For these parameters, we compare the ground-state phase diagram with that of 
the Anderson model and find a ferromagnetic region analogous to one found in the Anderson model, 
but which is pushed to somewhat higher densities and interaction strengths. In both models, we find 
a region within the ferromagnetic phase in which phase separation between a localized ferromagnetic 
domain and a weakly antiferromagnetic regime occurs. We then choose a set of parameter values 
appropriate for copper-oxide materials and explore the ground-state phase diagram as a function 
of the oxygen-oxygen hopping strength and the electron density. We find three disconnected re- 
gions of metallic ferromagnetism and give physical pictures of the three different mechanisms for 
ferromagnetism in these phases. 
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I. INTRODUCTION 



The nature of the microscopic description of metal- 
lic ferromagnetism, as found in transition metals such 
as iron, nickel and cobalt, is a long-standing problem 
in strongly correlated electron systems. The Hubbard 
model was formulated in the 1960'so in order to describe 
ferromagnetism in such materials. However, a ferromag- 
netic ground state in the single-band Hubbard model 
on nonfrustrated bipartite lattices has not been found 
at physical parameter values, and can be excluded in a 
large portion of the ground-state phase diagram. While 
mean-field theory yields a large region of stable ferromag- 
netism, fluctuations tend to destabilize the ferromagnetic 
phase, and numerical and variational calculations have 
narrowed the possible extent of a ferromagnetic state to 
a small region around the Nagoaka point. Therefore, it 
has become clear that additional features must be added 
to provide a description of metallic ferromagnetism.0a 

Recently, a number of such possible extensions to the 
Hubbard model have been investigated. In particular, 
there are three classes of additions which can enhance 
ferromagnetism: (i) a change in the noninteracting den- 
sity of states through the addition of frustrating hop- 
ping terms or the treatment of a geometrically frustrated 
lattice (ii) the inclusion of multiple orbitals per site or 
unit cell, and (Hi) the-addition of more general nearest- 
neighbor interactions na In relation to (i), Vollhardt et 
al.tj have emphasized that peaks or singularities in the 
occupied portion of the noninteracting density of states 
minimizes the kinetic energy loss due to polarization of 
the electrons. This picture is supported by numerical 



calculations in one@ and two dimensions,! and within 
the Dynamical Mean-Field Theory with a number of-dif- 
ferent forms of the noninteracting density of states B In 
addition, MiclkcO has shown rigorously that a Hubbard 
model with a less than half-filled flat band has a fully 
polarized ground state. This theorem has been extended 
to the case of nearly flat bands. For (ii), a number of 
multiband models, such as the periodic Anderson model, 
the Kondo lattice model, and the multiband Hubbard 
models with Hund's rule coupling, have been found to 
have ferromagnetic ground states.0 In relation to (Hi), 
Strack and VollhardttJ have studied a quite general form 
of the nearest-neighbor interaction and have found a fully 
polarized ground state at half-filling and one hole from 
half-filling (the Nagaoka state) under certain conditions. 

In this paper, we consider a one-dimensional model 
which contains elements of points (i) and (ii) from above. 
The model which we will study is a Cu-0 chain that 
contains one oxygen p-orbital and and copper d-orbital. 
Since one-dimensional models can usually be treated us- 
ing well-controlled analytical or numerical techniques, 
one can obtain a complete picture of the ground-state 
phase diagram without resorting to approximations nec- 
essary in two or three dimensions. For one-dimensional 
systems with nearest-neighbor hopping and an arbitrary, 
but real and particle-symmetric interaction, a theorem 
by Lieb and MattisO rules out a magnetized ground state. 
Because the Lieb-Mattis theorem no longer applies when 
electrons can pass around each other, the minimal addi- 
tion to the one-dimensional Hubbard model necessary to 
obtain ferromagnetism is a next-nearest-neighbor hop- 
ping. The resulting Hubbard model on a "zigzag" ladder 
has a ferromagnetic ground state in a substantial parame- 
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ter regimes One can also formulate more general, related 
models composed of triangular elements. Two such mod- 
els, studied by Tasakit2l and Penc et alJlil have been found 
to have ferromagnetic phases. The d—p model we study 
here is (aside from an added near-neighbor Coulomb re- 
pulsion) a particular case of the model studied in Rcf. 

Another approach to ferromagnetism in one- 
dimensional models, involving point (ii) from above, are 
models with multiple non-degenerate orbitals per site. 
The simplest of these models are the one-dimensional 
Periodic Anderson Model (PAM) and the Kondo Lattice 
model (KLM). Both of these models have ferromagnetic 



ground states in a relatively large parameter regime! 
These models can be related to one another in the regime 
of large /-repulsion for the PAM and small Kondo cou- 
pling for the KLM. The PAM is similar to the d-p model 
in that there are two linked non-degenerate levels per 
unit cell, and, as we will discuss in more detail in the 
following, exhibit similar behavior with the appropriate 
choice of parameters. 

Another motivation for studying the d-p model is the 
relationship with models for the CuC>2 planes of the high- 
T c superconductors. One dimensional d-p models have 
been thought to contain some of the essential ingredients 
needed to describe superconductivity in Cu02 planes.E-j 
Because of this, previous work on the one-dimensional d- 



p model has concentrated on determining the conditions 
under which superconducting cpsrelations are dominant 
for models witblia and without£3~ll3 direct hopping be- 
tween the p-orbitals. 

In this work, we will concentrate on the ferromagnetic 
phases in the ground state phase diagram and will discuss 
their physical origin. In Sec. [H], we introduce the models 
which we will stu dy an d discuss choices of the models pa- 
rameters. In Sec. Ill A , we compare the phase diagram of 
the Anderson lattice model from Ref. [l2| with that of the 
d-p model. We find that the ferromagnetic region per- 
sists in the d-p model but it is pushed to higher values 
of the interaction and higher fillings. There is a region 
of phase separation in which fer romag netic and antiferro- 
magnetic states coexist. In Sec. Ill B, we set the parame- 
ters of the d-p system to realistic values for the cuprates. 
We then investigate the phase diagram as a function of 
the oxygen-oxygen hopping amplitude, t pp , and the band 
filling, and find three disconnected ferromagnetic regions 
with three physically different mechanisms. 



II. MODELS 

We consider a one-dimension d—p system with Hamil- 
tonian 



J" 



£^£« - t P d £ (4tPj«t +Pl d ja) + Vpd £ nfr 



(1) 



<ij>a 



<tj> 



where denotes a sum over nearest-neighbor oxygen 
pairs, and (ij) a sum over copper-oxygen nearest neigh- 
bors. In the d-p system we work in the hole represen- 
tation so that d\ a {j>\ a ) creates a hole on d (p) site i 

(j) with spin a, and nf a — d itT d ia and n^ a are the local 
hole densities on a copper and oxygen site, respectively. 
The parameter t pp is the strength of the direct hopping 
between the p-orbitals, A is the difference in on-site ener- 
gies, U p and Ud are the on-site Coulomb repulsion on the 
p and d sites respectively, t p d is the hybridization and V p d 
the Coulomb repulsion between nearest-neighbor p and d 
orbitals. The lattice structure of the d p model including 
a schematic representation of the Hamiltonian parame- 
ters is shown in Fig. |l](a). Unless otherwise stated, we 
will consider lattices consisting of N copper sites and 
N + 1 oxygen sites, with open boundary conditions at 
the ends of the lattice so that the ends consist of oxygen 



sites with connections in only one direction. When the 
system has Nh total holes, we will discuss N p = Nh — N, 
which corresponds to the excess number of holes in the 
oxygen band when Ud > A > 0, and the corresponding 
p-band filling, n p = N p /N. 

This model has a gauge symmetry which can be used to 
permute the signs of the hopping matrix elements around 
each triangular element. For example, alternating the 
sign of t p d on succesive bonds and taking t pp positive is 
equivalent to having the same sign of t p d and taking t pp 
negative. Here we take the latter CcLSG, cLS shown in Eq. 
(0), restrict t p d to be positive, and allow the sign of t pp 
to vary. 

In order to understand the physics of ferromagnetism 
in the d-p model, it is useful to make a comparison with 
a model that has a similar structure in an appropriately 
chosen parameter regime, the one-dimensional periodic 
Anderson model. The PAM has Hamiltonian 

(2) 



H = - f £( c U+i<x + 4 +la c ia ) +£f j2nL + uJ2 4Ai + v £(4/- + ft*°) , 

la icr i ia 
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where cJ CT and Cj CT create and annihilate conduction elec- 
trons with spin a at lattice site i, and fj a and fi a create 
and annihilate local /-electrons. Here t is the hopping 
matrix element for conduction electrons between neigh- 
boring sites, £/ is the energy of the localized /-orbital, U 
is the on-site Coulomb repulsion of the /-electrons, and 
V is the on-site hybridization matrix element between 
electrons in the /-orbitals and the conduction band. We 
denote the number of electrons by N e , and N is the num- 
ber of unit cells (each consisting of one / site and one con- 
duction site) in the lattice. Since there are two electronic 
orbitals in each site, the quarter-filled case corresponds 
to N e — N and the half-filled case has N e = 2N. In anal- 
ogy to the d-p model, we will discuss N c = N e — N, the 
excess number of electrons in the conduction band when 
U > — £/ > 0, and the corresponding filling, n c = N c /N. 
Fig. 0(b) shows a schematic represention of the PAM. 

One can set the parameters of d-p model so that it is 
the same as the PAM except that the d-orbital hybridizes 
with the two nearest-neighbor copper atoms, while the 
PAM hybridizes only on-site. In order to compare the 
two models, we set V p d = U p = 0, t pp = t, t pi i = V, 
A = -e f and U d = U. 




U 

FIG. 1. A schematic diagram of (a) the copper-oxide lat- 
tice and (b) the periodic Anderson model, with the system 
parameters marked. 



III. CALCULATION AND RESULTS 

Since itinerant fcrromagnctism is an intrinsically 
strong-coupling phenomenon that occurs due to a suit, 
tie competition between kinetic and potential energy,EI 
it is generally quite hard to treat with analytic meth- 
ods. While exact statements can be Baade in some spe- 
cial cases, -such as the Nagaoka stateEJ and the case of 
flat bands Jj in general one must resort to mean-field the- 
ory or variational techniques, both of which give phase 
diagrams which can be at best qualitatively accurate. 
Therefore, sufficiently accurate numerical methods can 
be useful to determine the properties of models with 



ferromagnetic phases. Here we use the Density Matrix 
Renormalization Group to calculate the ground-state 
properties of the d-p model and the PAM. The DMRG 
is a variational numerical method closely related to ex- 
act diagonalizatkm. but which can be used to treat much 
larger systems .ElO Being able to treat sufficiently large 
systems is important since the models we study here 
have two fermionic sites, i.e. sixteen degrees of freedom 
per unit cell, which would severely limit the maximum 
size available for exact diagonalization. Quantum Monte 
Carlo methods would suffer from the fermion sign prob- 
lem at all band fillings on the d-p lattice. 

We use the finitepsystem version of the ground-state 
DMRG algorithmc2l to accurately calculate the energy 
expectation values of equal-time operators in the ground 
state including local spin, density, and various correla- 
tion functions. We keep up to 800 states in the system 
block and treat lattices of up to N — 32 unit cells. The 
maximum sum of discarded density matrix eigenvalues is 
approximately 5 x 10 . 

In order to search for ferromagnetic phases, it is essen- 
tial to be able to determine the total spin of the ground 
state. Since our version of the DMRG algorithm does 
not allow direct control of the total spin, we use a com- 
bination of methods to do this. First, we can calculate 
the expectation value (ip \S 2 \ipo) directly, where ipo is 
the variational DMRG ground state calculated in a par- 
ticular S z sector, and S 2 is the explicit operator for the 
square of the total spin. Since the DMRG ground state is 
variational, (ipo\S 2 \^o) does not always take on the quan- 
tized values due to mixing between states of different S z ■ 
Second, we can add a term of the form AS 2 to the Hamil- 
tonian, raising the eaergy of higher spin states within a 
particular S z sector. □ We can use this to calculate the 
ground-state energy as a function of S 2 by examining 
the energy of the minimum S* 2 ground state in a particu- 
lar S z sector. Finally, we can examine the ground-state 
energies as a function of S z and use the degeneracy in 
S z to determine S 2 . For application of these methods to 
the ferromagnetism in the PAM, see Ref. |lj and to the 
Hubbard chain with next-nearest-neighbor hopping, see 
Ref. |. 



A. Comparison with the periodic Anderson model 

In the first part of this work, we present a compari- 
son between the d-p system and the PAM. The quali- 
tative picture developed for the ferromagnetic phase of 
the PAMc3'll3 should also apply to the d-p model. In or- 
der to make the correspondence to the PAM, we choose a 
regime in which the (i-orbital lies below the bottom of the 
p-band. At exactly quarter filling, all the electrons are in 
the d-sites and all the spin states are degenerate if there 
is no hybridization. When the hybridization is turned 
on, the super-exchange favors antiferromagnetic interac- 
tions. The exchange coupling constant can be calculated 



3 



from perturbation theory giving 

1 

U d + 



J A = 



(A + v pd y 



2A + JX 



(3) 



or when U d = 2A and U p = V pd = 0, J A = 6^ d /A 3 . 
When a hole is added to the quarter-filled system, it will 
prefer to go to the p-sites since U d > |A|. This extra hole 



will tend to form singlets with the mostly localized spins 
in the d-sites. The hole can delocalize and thus lower its 
energy if the spins of the d-holes are oriented in the same 
direction, i.e. are ferromagnetically ordered. The binding 
energy for these singlets can be calculated within pertur- 
bation theory in a manner similar to that-applied to the 
two-dimensional case by Zhang and Rice,E3 yielding 



Js = 2f 



pd 



i 



l 



v pd 



V pd + U P Ud - A - V pd 



(4) 



For U d = 2A and U p = V pd = 0, J s = 8t 2 pd /A. 

In general, when there are N p — Nh — N holes in the 
p-band, this effect will favor a ferromagnetic ground state 
with total spin S = (N—N p )/2. We will denote this value 
of S for the ground state as complete ferromagnetism, 
while a value of S smaller but still greater than the min- 
imum, will be labeled incomplete. Note that a complete 
state has a lower magnetization than a saturated one. It 
represents a state where all the uncompensated d-spins 
are aligned. 

Thus, one would expect a competition between fer- 
romagnetism and antiferromagnetism near quarter fill- 
ing similar to that found in the PAM.L3 For the PAM 
there is an exact result that proves the ground state is 
ferromagnetic with complete, magnetization for N c = 1 
(N e =N+ 1) and U d = ooB@ There is no such result, 
however, for the d-p model. Furthermore, while the an- 
tiferromagnetic exchange is a sixth-order process in the 
Anderson lattice model, it is fourth order (and therefore 
stronger) in the d-p system. Therefore, it is not clear 
that the ferromagnetism will still be present in the d—p 
system. YanagisawaEa studied a 2 x 2 CuO cluster with 
exact diagonalization and found a ferromagnetic ground 
state. However, it is difficult to draw general conclusions 
about the behavior in the thermodynamic limit from such 
a small cluster. 

In our DMRG calculations, we consider chains with 
N = 16 and open boundary conditions and set t pp = 0.5, 
t pd = 0.375, A = U d /2 and U p = V pd = 0. This choice 
allows us to compare the phase diagram with the one of 
the Anderson lattice model from Ref. [l2] where t = 0.5, 
V = 0.375, e f = -U/2. In Fig. |, we present both 
phase diagrams in the plane of band filling and the local 
Coulomb repulsion, U or U d . Here n c = (n p = for 
the d-p model) represents the quarter-filled case, which 
in the large Coulomb repulsion regime means that the 
localized orbitals are singly occupied and the extended 
orbitals are empty in the zero hybridization limit. In the 
d-p system, Fig. |2j(b), the region of complete ferromag- 
netism is pushed towards higher values of U d and higher 
densities and there is a wide region of incomplete states 
at low densities. This is consistent with the fact that 
antifcrromagnetic correlations at quarter filling in this 
system are of lower order (and therefore stronger) than 
in the Anderson lattice model. 




(a) 



i ®©©© en i 




(b) 



FIG. 2. Ground-state phase diagram in the plane of band 
filling and the local Coulomb repulsion for an N = 16 sys- 
tem of (a) the Anderson lattice with V = 0.375, t = 0.5, 
€f = —U/2, and (b) the copper-oxide lattice with equivalent 



0.375, t pp = 0.5, A = U d /2, U p = V t 



pd 



o. 



parameters: t 

At n c = or n p = 0, each system is quarter-filled, i.e. has 
one hole per unit cell. The dotted and dashed lines in (b) 
correspond to Eq. (^) and Eq. respectively. 

In Fig. [3] we show the Cu-Cu correlation function, 
(Sf(r)S^(0)}, as a function of distance, r (measured in 
units of the lattice constant) , for the d-p model with the 
parameters described above and U d = 5 for N p = 0, 2 and 
4. Here S^(r) = d r1 .d r> i is the <i-spin raising operator 
on site r. For N p — (quarter filling), the correlations 
are clearly antiferromagnetic with the amplitude decay- 
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ing slowly with distance, while for N p = 2 and N p = 4, 
the correlations are ferromagnetic. This illustrates that 
the system goes from an antiferromagnetic state at quar- 
ter filling to a ferromagnetic state as soon as additional 
holes are added for these parameter values. 

The competition between the tendency to ferromag- 
netism and the antiferromagnetic exchange gives rise to 
the region of incomplete ferromagnetism at low doping. 
In this region, we find that the system phase separates 
into domains of complete ferromagnetism in which the 
p— holes are localized and antiferromagnetic domains in 
which the p-band is empty. Indeed, the first symptom of 
this phase separated state can be seen in Fig. |3|, in which 
there is a jump in the correlation function for r = 12 
for the case of 18 holes. This peculiar behavior is due to 
the tendency of the holes to localize and form a region 
of complete ferromagnetism near the center of the chain, 
leaving the p-orbitals empty near the ends. 



smaller than the lattice size. This behavior is also seen 
for Ud = 5 and Ud = 6. 
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FIG. 3. The Cu-Cu spin-spin correlation function as a 
function of distance, r, for the copper-oxide lattice with the 
same parameters as in Fig. |2] and Ud — 5 for 3 different fillings. 

To illustrate this phase-separated regime, we have cal- 
culated the chemical potential fx = E D (N/ l + 1) — E Q {Nh) 1 
where E (Nh) is the ground state energy with Nh holes, 
for chains of N = 8, 16 and 32 unit cells. We show \i as 
a function of the filling for Ud = 5 in Fig. [|. For 8 unit 
cells, (i increases monotonically with n p , but for 16 unit 
cells there is a shallow minimum at n p = 0.0625. For 32 
unit cells, there is a flat region which is shown in detail in 
the inset, with a minimum at n p = 0.09375 which is in- 
dicative of phase separation. We find a similar behavior 
of the chemical potential is found for Ud = 6 and Ud = 4. 

To illustrate the localization of the holes in the p-band, 
we display the density in the p-sites as a function of po- 
sition in the lattice for Ud = 4 in an N — 32 unit-cell 
chain with N p = 0,1,2,3 in Fig. ||. At quarter filling, 
the p-levels are almost empty. Upon doping, they begin 
to fill with holes, but instead of spreading out over the 
lattice, the additional holes localize in a region which is 



FIG. 4. The chemical potential n, as a function of the 
p-band filling, n„, for the Cu-O lattice with the same param- 
eters as in Fig. Ml for different system sizes. The inset is an 
expanded view of the low-density results for the 32 unit -cell 
lattice showing a minimum in the chemical potential, the sig- 
nature for phase separation. 

In Fig. ||we show the local spin density (Sf) as a func- 
tion of position in the lattice for an N — 32 chain with 
Ud = 5, N p = 2 and S z — S = 5 (the ground state) 
for both the p-band and the d-band. There are two dis- 
tinct domains in the lattice: one in which the (i-holes are 
completely polarized and another one in which they tend 
to order antiferromagnetically. The coupling between p 
and d holes is clearly antiferromagnetic. The asymme- 
try of the distribution is a numerical artifact due to the 
near-degeneracy of states with the phase-separated re- 
gion centered at different points. 

To gain some physical insight into the phase separation 
mechanism, one can use energetic arguments to estimate 
the stability of competing phases. The energy scale in 
the antiferromagnetic state is set by Ja, the antiferro- 
magnetic exchange energy, Eq. (||). The energy of the 
uniform ferromagnetic phase can be estimated by the 
binding energy of the Zhang Rice like singlets, JsN p . 
For a rough estimate of the magnetic energy scale, con- 
sider an antiferromagnetically ordered state versus one 
that has complete ferromagnetism. The Bethe Ansatz ep=. 
ergy of a one-dimensional Heisenberg antiferromagnetE3 
consisting of N — 1 bonds (measured relative to the ferro- 
magnetic state and neglecting end effects) is J A S (N — 1), 
where = ln2J^. In order for complete ferromag- 
netism to occur, the condensation energy of the singlets, 
JsN p must be larger than this, leading to the condition 

n p > In 2 Ja/Js = n P d ■ (5) 

This estimate is only valid for very low densities of p— 
electrons, since it neglects their kinetic energy, and for 
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strong coupling where the magnetic picture is valid. Ac- 
cording to this estimate, ferromagnetism is favored over 
antiferromagnetism for C/ d > 3 at all densities above 
quarter-filling on the lattice sizes that we have studied. 



However, it yields a smaller n pc \ than we find in the nu- 
merical calculations. For example, for Ud = 4 and the 
parameter set used in this section, n pc \ — 0.018. 
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FIG. 5. Density in the oxygen sites as a function of position in the lattice, i, with N = 32 unit cells, the same parameters as 
in Fig. | and U d = 4 for (a) N p = (quarter Ming), (b) N p = 1, (c) N p = 2, and (d) N p = 3. 



Now let us examine the possibility of phase separation 
between a region of undoped antiferromagnetism and a 
region of complete ferromagnetism. In order for such a 
phase to exist, the energy gained by forming an undoped 
antiferromagnetic domain must compensate the loss in 
kinetic energy of the Zhang-Rice singlets, which have to 
localize in a region smaller than the lattice size N. We 
estimate the kinetic energy of the Zhang-Rice singlets 
by taking them to be noninteracting hard-core bosons 
confined to a lattice of length X with open boundary 
conditions. (We study lattices with open boundary con- 
ditions here, so that this form is valid up to the number 
of p-sites, X = N + 1.) The single-particle energy levels 
are given by 



Ej = —2tzRCOS 



X + l 



(6) 



where tzR is the effective hopping of the singlets and 
j = l,...,N. Since hard-core bosons are equivalent 
to spinless fermions in one dimension, the total kinetic 
energy is given by successively occupying single-particle 
states. 



T(N p ,X) = J2^ 



dj Ej 



X 



N p TT 



-2t Z R — sin 

7T A 



(J) 



G 



where the approximation on the right-hand side is valid 
in the continuum limit: X ^> 1. The energy gained by 
constricting the ferromagnetic domain to a length X is 

E PS (N p ,X)= jf(N-l-X) 

+ T(N p ,N + l)-T(N p ,X) . (8) 

The phase-separated state will be favored over the ho- 
mogeneous ferromagnetic state if Eps is positive. In that 
case, the size of the ferromagnetic region is given by the 
value of X that maximizes the energy gain, determined 
by 

- J a + 3X 3 . ( 9 ) 

where we have taken Npir <C 1 in order to expand the 
sine and have used the continuum expression from the 
right-hand side of Eq. . Eq. (^) has one maximum for 
X > which occurs at 



* max = n p (jjjar) = n p x o • ( 10 ) 

It is easy to verify that Eps(N p , X max ) is always posi- 
tive as long as < X max < N. Therefore X max is the 
length of the ferromagnetic domain when phase separa- 
tion is present. Note that X max oc N p , in agreement with 
the behavior seen in Figs. ^| and |6|. In order to numeri- 
cally estimate Xq for the d-p system, we take tzp = t pp 
and J^ ff = ]h2Ja, where J a is given by Eq. (||) for the 
d-p model. (One could calculate tzn more accurately in 
perturbation theory, but only its order of magnitude is 
important, since it appears within the cube root.) For 
the parameters used in this section, this gives xq — 7, 8.5 
and 10 for Ud = 4, 5 and 6, respectively. This is in rea- 
sonable agreement with our numerical results (See Fig. 
^) , which also show an increase in the localization length 
with Ud- 
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FIG. 6. The local spin density (Si) as a function of lattice position, i, for (a) the oxygen sites and (b) the copper sites on 
an N = 32 lattice with S z = S = 5 (the ground state) and N p — 2. Here Ud = 5 and the remaining parameters are the same 
as in Fig. 



The boundary between the uniform ferromagnetic 
phase and the phase-separated phase is determined by 
the condition X max = N , leading to a critical density 

( 3J A 

The densities n pc \ and n pC 2 calculated from Eqs. (j^) and 
( pd| ) are represented in the phase diagram of Fig. ||(b) 
as dotted and dashed lines, respectively. Notice that the 
numerical calculation yields a somewhat wider region of 
phase separation than given by these estimates. In order 
to estimate the finite-size effects, we have also examined 
the spin of the ground state on N = 8 and N = 32 unit- 
cell chains. We find that the points of complete polar- 
ization in Fig. ||(b) behave consistently with system size, 



i.e. are also completely polarized for N = 8 and N = 32. 
However, the scaling behavior of the incomplete states 
is more complicated. In particular, some states which 
are incompletely polarized in the N — 16 chain become 
completely polarized for N = 32. These cases are circled 
with a dashed line in Fig. ||(b). Therefore, the phase- 
separated region tends to become narrower in the ther- 
modynamic limit. There is, however, evidence that some 
of the incompletely polarized states remain in the ther- 
modynamic limit. For Ud — 5 and 6, the two densities 
closest to quarter filling on the N = 16 chain have total 
spin S — 5/2 and 5, and S = 7/2 and 6, respectively. On 
the N = 32 chain, the ground-state S at the same den- 
sity doubles, suggesting that S/N will remain constant in 
the thermodynamic limit. At Ud — 4 and Ud = 3, none of 
the incomplete states show such a simple scaling and we 
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believe that larger system sizes are needed to determine 
the nature of the phase boundaries. For Ud = 3, we find 
no evidence of phase separation at any finite system size. 
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FIG. 7. (a) The chemical potential, versus density, n c , 
for the periodic Anderson model with V = 0.75, t — 0.5, and 
€f = — (7/2 on an N = 32 lattice. It shows a minimum for 
all three [/-values, (b) Conduction electron density versus 
site for the (7 = 5 case, showing localization of the excess 
conduction electrons in a region smaller than the lattice size. 

Let us now examine the relevance of the phase sepa- 
ration arguments to the PAM. Our energetic arguments 
should also apply with the appropriate mapping of J a 
and Js to the model parameters. For the symmetric 
case, e s = -U/2, J s = 8V 2 /U and J A = %VH 2 /e). For 
the parameter sets we have compared on the two models, 
J A is smaller in the PAM than in the d-p system. This 
implies that phase separation occurs over much larger 
localization lengths in the PAM. Indeed, for the conduc- 
tion electron density n c = 1/32 we obtain x a » 28, 21 
and 14 for U — 6, 5 and 4, respectively. Because these 
lengths are larger than those for the d — p system, we 
could not find definite numerical evidence of a phase- 
separated state on the lattice sizes that we were able to 
study; larger systems sizes would be required. In order to 
confirm our picture of the origin of the phase separation, 



we can instead reduce Xq by increasing the hybridization 
V; we double V to V — 0.75, decreasing the localiza- 
tion length by a factor of two to x w 14, 10 and 7 for 
U = 6, 5 and 4, respectively. Using the DMRG, we then 
examine N — 32 systems with t — 0.5, e/ = —U/2 and 
V = 0.75 for U = 4, 5, and 6. The results, presented 
in Fig 0, show numerical evidence for phase separation 
for all three [/-values. The chemical potential, Fig. 0(a), 
has a minimum as a function of the filling, and the con- 
duction electron density, Fig. 0(b), tends to localize in 
a region smaller than the size of the system. Therefore, 
we find that phase separation can also occur in the PAM 
under the right conditions. 
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FIG. 8. Phase diagram in the plane of oxygen occupation, 
N p , and hopping t pp for an N — 16 copper-oxide lattice with 
t p d = 1, U d = 8, A = 3, Up = 4, V pd = 0.5. The total spin, 
S, has been calculated everywhere in the plane, but is shown 
as a numerical value only when it is greater than the mimi- 
mum value, 5* = for even or S — 1/2 for odd N p . Cases of 
"complete" polarization are circled. 



B. Phase Diagram for "Realistic" Parameters 

In the second part of our work, we think of the d— 
p system as the one-dimensional analog of the copper- 
oxide planes in high temperature superconductors. We 



set the parameters t p 



1, U d = 8, U p 



,= 3 



and Vpd = 0.5, values reasonable for the cupratesJiH and 
then vary t pp and the hole filling. The resulting phase 
diagram is shown in Fig. ^ for an N = 16 chain and 
both possible signs of t pp . The cases with complete po- 
larization, S — (N — N p )/2, are circled. For t pp > the 
phase diagram has one ferromagnetic region consisting 
almost exclusively of complete polarization for t pp > 1 
and from n p ~ 0.4 to n p ~ 0.75. The arguments of the 
first part of our work should also apply to this region, 
i.e. the ferromagnetic ordering of the /-spins is caused 
by the condensation of Zhang-Rice singlets. As can be 
seen from the lack of regions of incomplete polarization, 
we find no phase separation at these parameter values. It 
is clear that Eq. (||) cannot properly predict the bound- 
ary of the ferromagnetic phase since it is invalid at the 
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relatively large N p values at which we find the phase. 
From Eq. ( |l0| ) , the charactistic size of the localization of 

1/3 

the singlets is given by xq w 4.5 t pp for these parame- 
ter values. Our energetic arguments therefore predict no 
phase separation at the densities at which we do find a 
ferromagnetic phase. 

When tpp < the phase diagram is quite different. 
There is one point at quarter filling and t pp — —2 for 
which the ground state is fully polarized (S = 8 for 
N = 16), surrounded by a region of smaller polariza- 
tion. Additionally, there is a region with polarization 



greater than complete, but less than fully polarized for 
a narrow range of t pp but a wide range of hole occu- 
pation, N p . More precisely, the extent of this region is 
-0.5 < tpp < -0.7 and 7 < N p < 14. 

We first discuss the fully polarized phase at quarter fill- 
ing and tpp ks —2.0. One can gain insight into the origin 
of this ferromagnetic phase by examining the noninter- 
acting band structure. Diagonalization of the Hamilto- 
nian, Eq. ([l]), with Ud = Up = Vpd = and periodic 
boundary conditions leads to the hybridized bands 



(12) 



e ± (Jb) = A/2 



tpp cos k I 



^(A/2 - tpp cos fc) 2 + 2f + cos k) 



These bands, plotted for three t pp values near the fully 
polarized point, are shown in Fig. ^. As can be seen, the 
lower band does not overlap with the upper band and has 
a dispersion that is quite flat; it becomes completely flat 
at tp P — —1.78. Therefore, all the holes will go into ttue, 
half-filled flat band at this point. A theorem by Mielke,u 
subsequently extended to nearly flat bands by Mielke and 
Tasaksfl shows rigouously that for the Hubbard model 
with a low-lying completely flat band at half-filling or 
less will have a fully polarized ground state. Therefore, 
this region is due to flat-band ferromagnetism. This is 
substantiated by the strong dependence of the polariza- 
tion on tpp and n p ; there is only a small region of partial 
polarization surrounding the fully polarized point. 




FIG. 9. The dispersion of the non-interacting bands, 
e±(fe), of the d-p lattice for t p d — 1.0, A = 3.0 and three 
different t pp . 

We next treat the narrow strip of ferromagnetism at 
tpp « —0.5 and intermediate N p . The first important 
issue is whether the partial polarization in this region 



persists in the thermodynamic limit. In Fig. [T(J, we plot 
S/N as a function of 1/N for three different fillings. As 
can be seen, S/N increase linearly with 1/N and extrap- 
olates to a finite value as N — > oo. 
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FIG. 10. Total spin per site, S/N, as a function of 1/N for 
a copper-oxide lattice with the same parameters as in Fig. 
[s[ t pp = — 0.5 and three different fillings. The corresponding 
solid symbols show the N — > oo values obtained from a linear 
extrapolation in 1/N. 

One can gain some insight into the relationship be- 
tween the sign of t pp and the mechanism for ferromag- 
netism by examining the behavior of small clusters. In 
particular, we examine the behavior of one triangular ele- 
ment of the lattice, i.e. a cluster consisting of one copper 
and two oxygen sites, occupied by two holes (N = 1, 
N p = 1). This Cu02 cluster can be treated analytically 
if spin and reflection symmetry are taken into account: 
The total spin can either be zero or one, and states can be 
symmetric and antisymmetric with respect to exchange 
of the oxygens. The antisymmetric 5 = and the sym- 
metric S = sectors can be represented by 2 x 2 matrices, 
the symmetric S = sector is 4 x 4, and the antisymmet- 
ric S = 1 sector is 1 x 1. All of these matrices can be di- 
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agonalized analytically (although it is easier, in practice, 
to treat the 4x4 matrix numerically since the general 
solution is rather complicated). Taking t p( i = 1, U4 = 8, 
Up = 4, A = 3 and V p d — 0.5, and varying t pp , we find 
that for tpp > 0, the ground state is always a singlet 
with positive reflection symmetry. For t pp < 0, there is 
a competition between singlet formation, favored by di- 
rect antiferromagnetic exchange, and triplet formation, 
favored by third order exchange around the plaquette, 
which is ferromagnetic. 
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FIG. 11. The energies of the three lowest lying states of 
the C11O2 cluster for t p d = 1, Ud = 8, U p = 4, A = 3 and 
Vpd — 0.5, as a function of — t vv . Here "SYM" denotes a state 
that is symmetric with respect to exchange of the oxygens, 
and "ASY" denotes an antisymmetric state. 



point of maximum polarization occurs when n p ~ 0.5, i.e. 
there is one hole for every second p-site. At this filling, 
ferromagnetism based on exchange around the triangular 
elements would be somewhat frustrated since each p-site 
is shared by two triangular elements. These mechanisms 
would tend to suppress ferromagnetism. Also, as seen in 
Fig. [ll], the triplet ground state is close in energy to a 
singlet excited state even on a single triangular element. 
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FIG. 12. Local spin density (Sf), as a function of Cu-site 
index i for both the copper (Cu) and oxygen (O) sites on an 
TV = 16 d-p lattice with N p = 8 (n p = 0.5), t p „ = -0.6, 
Sz = 5 and the remaining parameters as in Fig. pi For the 
oxygen sites, i takes on half-integer values between the corre- 
sponding Cu-sites. 



In Fig. [ll], we show the dependence of the energy of 
the three lowest lying states of the CuC>2 cluster on —t pp . 
As can be seen, for < — t pp < 0.46 and — t pp > 4.5, the 
ground state is a singlet with positive reflection symme- 
try. In the intermediate range, 0.46 < — t pp < 4.5, the 
ground state is an antisymmetric triplet, which is close 
in energy to the first excited state, an antisymmetric sin- 
glet. This illustrates the effect of the sign of t pp on the 
basic triangular element of the lattice: a localized triplet 
is only possible for negative t pp and can occur only for 
a certain range of — t pp , when there are two holes on the 
plaquette. 

While the triplet ground state of the CuC>2 cluster mo- 
tivates the possibility of a ferromagnetic state in this re- 
gion, and provides an explanation for the minimum value 
of — tpp required, it cannot provide a complete picture of 
this phase. In particular, the range of —t pp with a triplet 
ground state is much wider in the CuC>2 cluster than 
the range of —t pp in the ferromagnetic strip in Fig. ||. 
When more triangles are connected together, interaction 
between holes in the p-band can become important. The 



We can gain additional insight into these effects from 
diagonalization on clusters larger than three sites. For 
clusters of N = 2, 3, and 4 unit cells (5, 7, and 9 to- 
tal sites) with open boundary conditions, a ferromag- 
netic ground state appears at all fillings in the range 
< N p < N for -tpp > -<™ n , with i™ n ranging from 
-0.52 to -0.48. However, the extent of this region in t pp 
and the value of the polarization depends on N p . For 
N p = 1, one hole in the p-band, and all three sizes, we 
obtain a ground-state spin at or close to the maximum 
possible value, up to £^f x ranging between -2.9 and -4.2. 
This behavior is like that in the N = 1 CuC>2 cluster, 
since all the d-holes are ferromagnetically polarized. It 
is also reminiscent of the ferromagnet phase in the PAM, 
except that it is a triplet state which delocalizes rather 
than a singlet. Similar ferromagnetic states have been 
found in the one-dimensional Kondo lattice model with 
ferromagnetic exchanged 

For larger N p , we obtain a ferromagnetic ground state 
up to £™ ax ranging from -0.62 to -0.70, a value in good 
agreement with that found in larger systems. The 



10 



ground-state spin varies and depends on N p . One can un- 
derstand the nature of this phase by examining the distri- 
bution of spin polarization and the spin-spin correlation 
function. We display the local spin polarization, (Sf), 
for parameters that yield the maximum polarization on 
a N = 16 cluster in Fig. [l2|, taking S z — 5, the maxi- 
mal value. One can see that the holes on the d-sites are 
antiferromagnetically aligned, which we have confirmed 
by examining the spin-spin correlation function. The 
p-sites, however, all have positive polarization, and are 
ferromagnetically correlated with the positively polarized 
d-sites, i.e. with every other d-site. The correlations be- 
tween p-sites are generally ferromagnetic. The picture 
is then that there is an alternation between triangular 
plaquettes with a ferromagnetic moment and isolated d— 
holes antiferromagnetically aligned with the neighboring 
plaquettes. This phase is reminiscent of a ferrimagnet 
with alternating S — 1/2 and S — 1. Changing the band 
filling can frustrate this order, changing the amount of 
polarization. One domain wall is visible in Fig. 12. Since 



this phase is due to a subtle competition between ferro- 
magnetic and antiferromagnetic exchange processes, it is 
not surprising that it is quite sensitive to the value of t pp . 
Certainly the picture of localized plaquettes breaks down 
as the hopping between the p-sites increases. 



IV. CONCLUSIONS 

We have studied ferromagnetism in the phase diagram 
of the one dimensional d—p model. In the first part of 
this work, we have compared the phase diagram near 
quarter filling of a d-p system to the one-dimensional pe- 
riodic Anderson model. In both cases, we have found a 
ferromagnetic region near quarter filling, but, for equiv- 
alent values of the parameters, the ferromagnetic region 
in the d-p system is pushed towards higher values of the 
Coulomb interaction and higher densities. 

Between the antiferromagnetic state at quarter filling 
and the homogeneous ferromagnetic region, there is a 
phase-separated regime in which ferromagnetic domains 
are segregated from the antiferromagnetic background. 
The phase-separated state arises because the loss of ki- 
netic energy of the holes due to their localization in a 
small region is compensated by the gain in antiferromag- 
netic bonds. This regime is more evident in the d-p sys- 
tem because the antiferromagnetic coupling is stronger 
than in the Anderson model and therefore the charac- 
teristic length for the ferromagnetic domains is smaller. 
Nevertheless, we have presented numerical evidence for 
the phase-separated state for both models. 

Recently, similar phases with a hole-rich ferromag- 
netic domains segregating from hole-undoped antiferro- 
magnetic regions have been reported for models such as 
the one-orbital ferromagnetic Kondo model and other 
models studied in relation with manganites. (See Ref. 
|28| and references therein.) In our case, it must be em- 



phasized that the ferromagnetism is not saturated as in 
the case of these models. In the case of the ferromag- 
netic Kondo model, the phase-separated state was found 
near half-filling. However, it was pointed out that such a 
phase-separated state can arise at low densities if an anti- 
ferromagnetic Heisenberg coupling is added to the model. 
Several studies in different numbers of dimensions using 
different methods suggest that the phase-separated state 
is very robust. It is also mentioned that the inclusion 
of longer-range Coulomb repulsion can turn the phase- 
separated state into more exotic phases, such as striped 
phases with ferromagnetic domains. While we have only 
studied a one-dimensional system, the arguments pre- 
sented for the physical mechanism of phase separation 
should still be valid in higher dimensionsal systems. 

Also, it is important to notice that the Anderson model 
can be mapped to the antiferromagnetic Kondo model for 
large values of the Coulomb repulsion. It is well known 
that a ferromagnetic phase is also present in the an- 
tiferromagnetic Kondo lattice model.E3 However, there 
is no phase-separated state because the antiferromag- 
netic coupling is higher order than the mapping. We 
believe that an additional antiferromagnetic Heisenberg 
coupling of appropriate strength between the localized 
spins would induce a phase-separated state in the Kondo 
lattice model. 

In the second part of our work, we have studied the d-p 
model the parameters to values relevant to the cuprates. 
We have found that the sign of the oxygen hopping t pp 
is crucial in determining the physics of the ground state. 
Three disconnected ferromagnetic phases were found in 
the phase diagram. One phase occurs at positive t pp and 
is similar to the ferromagnetic phase found in the PAM 
and discussed in the first part of this work. A second 
phase appears for negative t pp at and near quarter fill- 
ing and occurs at a point at which the occupied nonin- 
teracting band has a flat dispersion. This "flat-band" 
mechanism ia_i» agreement with a picture based on ex- 
act theoremsfl'Q A third ferromagnetic phase appears for 
small negative t pp and it is related to ferromagnetic elec- 
tron exchange processes on plaquettes with a positive 
product of hopping integrals. This phase has elements of 
fcrrimagnctism due to an ordering of triangular plaque- 
ttes with a ferromagnetic moment and isolated, antifer- 
romagnetically aligned d-holes. 
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